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Abstract The subsurface structure of an “average” supergranule is derived from 
existing HMI pipeline time-distanc e data products and compar ed to the best 
helioseismic flow model detailed in Duvall and Hanasoeel ( 20131 ). We find that 
significant differences exist between them. Unlike the shallow structure predicted 
by the model, the average HMI supergranule is very extended in depth, exhibiting 
horizontal outflow down to 7-10 Mm, followed by a weak inflow reaching a depth 
of ~ 20 Mm below the photosphere. The maximal velocities in the horizontal 
direction for the average supergranule are much smaller than the model, and 
its near-surface flow field RMS value is about an order of magnitude smaller 
than the often-quoted values of ^ 250 — 350 ms“^ for supergranulation. Much 
of the overall HMI supergranule structure and its weak flow amplitudes can be 
explained by examining the HMI pipeline averaging kernels for the near-surface 
inversions, which are found to be very broad in depth, and nearly identical to one 
another in terms of sensitivity along the z-direction. We also show that forward- 
modeled travel times in the Born approximation using the model (derived from a 
ray theory approach) are inconsistent with measured travel times for an average 
supergranule at any distance. Our findings suggest systematic inaccuracies in 
the typic al techniques used to s tudy supergranulation, confirming some of the 
results in Duvall and Hanasoed ( 2013l l. 
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1. Introduction 


Supergranulation has been an active area o f sola r research now for several 
decades. Since its initial discovery by iHart ( 1954 1. extensi ve work ha s been 
done to char acterize these structures at the photosphere le.g.lHart.ll956l: Simon 
and Leighton. II964 Worden and Simonl. 19761: Hathaway et al . 2oCT 2002 : 


Hirzberger et ajJ . 20f)8l : lRientord et all 2008 ; Roudier et al . , 2014 ; Williams et al.. . 

20141 . among others!. However, many open questions still remain regarding their 
subsurface flow structure and vertical extent through the convection zone. To 
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better characterize the subsurface properties of supergranulation, the methods 
of local helioseismology must be employed. 

One of the fi rst studies of supergranulation from a helioseismic standpoint was 
carried out by iKosovichev and Duvalll ( 1997f ) , where large-scale flow structures 
were observed through time-distanc e invers i ons do wn to a depth of 5 Mm below 
the photosphere. Follow-up work by Duvalll ( 1998ll suggested that supergranules 
exhibit a radial outflow down to a depth of ~ 8 Mm, followed by a deeper inflow. 
These findings were based on two-dimensional spatial correlations of near-surface 
inversions with flowsrecoveredat^|rious_de 2 thSj_Subseauentjnversions^^Zhao 
and Kosovichev (1200.81 1 also r eported the presence of an infl o w at depths of 9- 
12 Mm. A later analysis by Ijackiewicz. Gizon. and Birch ( 2008 ) found that 
supergranulation persists down to at least ^ 5 Mm, but it was determined that 
the signal at larger depths is completely dominated by random noise and is 
therefore i rretri evable with short data sets of duration ~ 24 hrs. Additionally, 
Woodard ( 20071) found that inversions below depths of 4 Mm were infeasible due 


to high levels of noise. 

To quantify the accuracy of these result s, some effo r ts hav e been made to 


validate current time-distance techniques. Zhao et al. ( 20071) performed the 


first time-distance analysisofjM 2 e;Scale,^^merica^,si^mlatio^,dat^^o^Ben- 
son, Stein, and Nordlund ( 200fi1 ). It was determined that it was not possible 
to accurately recover horizontal flows below a depth of ~ 3 Mm, and near¬ 
surface inversions for the vertical flow component were actually anticorrelated 
with those of the simulationjJikel^,due,^o_crossAa]k£ffiects^,^ieceny^DeGrave, 
Jackiewicz, and Rempel ( 20141) performed a full time-distance analysis of two 
1 day magne tohydrodynamic quiet-Sun simulations (for simulation details, see 


Rempell . 2014h . each containing large-scale supergranule-sized flows, and found 


that recovered flow coherence was lost at depths exceeding 3-5 Mm. In fact, 
inversions performed at depths of 7-9 Mm actually suggested the presence of 
spurious flows that we r e not present in the si mulation. This result has since been 
suppo rted by Svandal (201^. Inversions bv IPeCrave. Jackiewicz. and Remr^ 
( 20141 ) for the vertical flow component were also found to be unreliable at all 
depths, even after trying to minimize the effects of cross-talk. 

It is becoming more apparent that recovering the subsurface flow properties 
of individual supergranules through helioseismic inversions at depths greater 
than 5 Mm is quite difficult due to high levels of random noise. To this end, 
there have been several attempts to study these features from a statistical 
viewpoint. Averages can be ta ken over many th ousands of features to help 
boost the signal-to-noise ratio ( Birch et all . l20()6[) . To date, these analyses of 
an “average” supergranule have relied on identifying regions of strong out¬ 
flow located in near-surface divergence (or nearly equivalently, point-to-annulus 
travel-time difference) maps. These near-surface av erages have then been u sed 


to more accurately measure surface properties (e.g. Duvall and Birchl. 20ld) or 


to constrain subsurface flow models ( Duvall and Hanasogel . 20131 : Svanda . 2012 


Duvall. Hanasoge. and Chakrabort 


Recently, iDuvall and Hanasogi 


i (0 


20141) of supergranulation. 

2013h derived a series of supergranule flow 
models based on averaging features measured using large-distance separation 
(A = 2-24°) point-to-annulus travel-time difference maps. It was observed that 
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Figure 1. A cut in depth at y = 0 through the DH 2 supergranular flow field, where arrows 
denote the vectors flow in the Vx — Vy plane. The longest arrow represents about 700 nis“^. 
The contour lines mark the 100, 200, 400, and 600 ms“^ scalar velocity |u| levels corresponding 
to the color bar. 


the measured travel-time differences at the center of the average supergranule 
approached a nearly-constant 5.1 sec for the largest A values. The idea is that 
waves propagating over these large distances penetrate so deeply through the 
convection zone that their paths are nearly vertical when leaving and approach¬ 
ing the solar surface, and are therefore likely to be mostly sensitive to the v ertical 
component of the supergranular flow field. Diivall and Hanasogd ( 201,ll l argue 
that the current value for the s urface vertical velocity component of supergran - 


ulation, = 0) ~ 10 ms ^ ( Duvall and Birch . 2010l : Hathaway et al. . 2002 1 


is too small to account for this 5.1 sec travel-time measurement. Therefore, the 
vertical flow field must increase in magnitude and peak below the surface at 
some depth. Constrained by the 10 ms“^ surface Vz measurement and assuming 
a Gaussian profile in the vertical direction, along with a mass-conserving flow 
held, they hnd the model that best approximates the data (referred to as DH 2 
hereafter) has the following properties: 


DH2: 

Vz{z = 0) = 10 ms-i and = -2.3 Mm) = 240 ms-\ FWHM^ = 2.1 Mm 

^;h(z = 0) = 138 ms-1 and v^^^{z = -1.62 Mm) = 700 ms"! 


for the vertical (vz) and horizontal (wh) how components, respectively. Figure[I] 
shows a cut through the DH 2 supergranule {vx, Vz) how component. Such a shal¬ 
low feature might appear quite surprising in light of the various works discussed 
previously. 

The goal of this work is to test the Duvall and Hanasoee ( 2013h how model 
against an “average” supergranule derived using existing helioseismic data prod¬ 
ucts provided by the Joint Science Operations Center (JSOC) Data Resource 
Management System (DRMS) for the HMI time-distance analysis pipel i ne (fo r 


details about these data products, see Couvidat et a7l . [ 2 OI 2 I : IZhao et al. . 2012h . 


This will allow us to directly examine the subsurface properties of an average 
supergranule at depth to compare to the results from these novel measurements. 
In Section[2l we briehy describe the HMI data and their properties. In Section[3l 
we explain the algorithm employed to identify and average supergranules in these 
data. Section!?] describes the resulting HMI supergranule in terms of travel-time 
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Table 1. Short-hand notation for the various flow 
features described throughout this paper. 


feature 

abbreviation 

HMI average supergranule 

HMI 

D&H (2012) Model 2 

DH 2 

Rempel average supergranule 

RSGi 

Rempel Model 

RSG2 

Best-fit Model (short-distance) 

BFMi 

Best-fit Model (large-distance) 

BFM2 


Table 2. Parameters for each of the eight HMI phase-speed 
data filters (adapted from IZhao et aZH2012D . 


filter No. 

annulus range 
[heliographic degree] 

phase-speed 

[kms“^] 

FWHM 

[kms“l] 

1 

0.54-0.78 

14.87 

4.37 

2 

0.78-1.02 

17.49 

4.37 

3 

1.08-1.32 

21.43 

5.47 

4 

1.44-1.80 

28.83 

9.40 

5 

1.92-2.40 

36.48 

5.91 

6 

2.40-2.88 

40.62 

5.17 

7 

3.12-3.84 

47.33 

8.29 

8 

4.08-4.80 

55.94 

8.95 


and flow fleld averages. Forward modeling comparisons between various super¬ 
granule models in the Born approximation are made in Section [5] Concluding 
remarks are made in Section [71 and a short discussion of supergranule mass 
conservation is given in Appendix |21 To avoid confusion, we note that the short¬ 
hand notation for the various flow features used throughout this work is given 
in Tabled] 


2. Data 

To briefly summarize, the HMI pipeline divides each raw, full-disk, 8 hr Doppler 
velocity time series into twenty-five smaller patches, each spanning 30° x 30° 
and centered at 0°, ±24°, and ±48° in solar latitude and longitude. These data 
are tracked and remapped using Postel’s projection, and subsequently filtered 
in Fourier space using a series of Gaussian phase-speed filters (referred to as 
tdi-tdg hereafter, parameters of which are given in Table ID). Cross-covariances 
are measured from the filtered data in the point-to-annulus and quadrant con¬ 
figurations over a range of annuli radii (A). Travel-time differences are then 
computed from the cross-covariances in the out-in (‘oi’, denoted i5toi), west-east 
(‘we’, denoted (5twe), and north-south (‘ns’, denoted dtns) geometries (Duvall 
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et al. 

. I 997 II using two definitions: a “least-sauares” method defined bv Gizon 


DuvalL 199^ Couvidat_e|ja^ 2012h . For the current work, we are only interested 
in the Gizon and Birchl ^ 2002 1 ‘ oi’ measurements. Each travel-time map is of size 


256 X 256 pixels, with horizontal grid spacing da; = dy = 0.12° pixel”^. Travel 
times are then inverted for the horizontal flow components {vx-,Vy) in the Born 
approximation, centered around eight depths (0-1, 1-3, 3-5, 5-7, 7-10, 10-13, 
13-17, 17-21 Mm). Travel-time maps computed using several additional filters 
of higher phase-speed, as well as inversion flow maps centered at several depths 
larger than the eight mentioned here are also available from JSOC. However, 
these have been deemed unreliable by the HMl pipeline team, and have not 
been fully tested. They are therefore not included in the subsequent analysis. 

JSOC allows one to download travel-time and velocity maps measured within 
some user-defined time frame, centered at any of the twenty-five desired latitude, 
longitude pairs. Here, data were selected to cover the 32 days spanning 10 June 
2010 - 11 July 2010 (i.e. the same time span analyzed bv iDuvall and Hanasoge . 
I 2 OI 3 I I. Only the nine patches centered at 0°, ±24° were used to avoid the effects 


of foreshortening that one might encounter toward the limb ( Zhao et ai . 2012h . 


These areas overlap, some as many as four times. We also use only data computed 
from the first 8 hr time series of each date (00:00-07:59 UT). Since supergranules 
only live for ~ 1.5 days on average, individual cells will have evolved significantly 
from one time series to the next, and we can effectively treat these as different 
features. In all, we are left with a total of 72 travel-time maps per date (one map 
for each of the eight filters over the nine 30° x 30° patches), as well as 72 flow 
maps for each of the three flow components. 

To aid in the analysis, we also have available the line-of-sight (LOS) magne¬ 
tograms covering this same 32 day period. These allow us to identify any areas 
of strong magnetic field that might be present in our data to avoid averaging 
supergranules near these regions. 


3. Supergranule Identification 


An algorithm was written to identify supergranules and subsequently average 
them about all cell-center locations. The supergranules are identihed directly 
from the ‘oi’ travel-time maps. The algorithm begins by reading in the 72 travel¬ 
time maps for a given date. These “raw” maps are in the Hostel projection. To 
properly overlay the maps and average their overlapping regions, all travel times 
were first projected back into latitude and longitude coordinates. From here, 
maps corresponding to a particular phase-speed filter were then combined via 
interpolation onto a large Cartesian grid of spacing 1.4577 Mm pixeP^. This 
results in a total of eight large travel-time maps (one for each phase-speed filter) 
for a particu lar date, an example of which is shown in Figure [2j 

Following iDnvall and Hanasogel ( 2013h . the large tdi (lowest phase-speed) 
map was smoothed via convolution with a two-dimensional Gaussian function of 
(7 = 4 Mm, where the Gaussian FWHM = 2(7-\/2ln2. This serves to smooth out 
most of the small-scale noise fluctuations while preserving the overall large-scale 
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Figure 2. An example tdi ‘oi’ travel-time map for the 10 June 2010 date. This results from 
combining the nine 30° x 30° patches centered at 0°, ±24° and remapping onto a Cartesian 
grid. The color scale has been clipped to more easily see the regions of presumed large-scale 
outflow (negative ‘oi’ travel time) which trace the positions of supergranule cells. 



X [Mm] X [Mm] 


Figure 3. Left: A close-up of the central region of Figure[2](in negative) after smoothing via 
convolution with a 2D Gaussian function of cr = 4 Mm. Supergranule cell centers identified by 
our algorithm are marked with black asterisks. Right: The cell centers (now marked by white 
asterisks) plotted over the same area of the unsmoothed travel-time map. 


supergranule structures. A simple algorithm was developed that takes each tdi 
travel-time map and finds the coordinates of large negative signals vertically 
(column-wise) and horizontally (row-wise). It then determines the coordinates 
common to both searches, and discards candidate cells separated from one an¬ 
other by less than 23 Mm. Any candidate within 23 Mm of the domain edges 
is also removed to avoid the counting of partial cells. Lastly, regions of strong 
magnetic field (-Blos > 250 G) were located using the magnetograms, and any 
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cells within these areas were eliminated. As an example, Figure [3] (left) shows 
the central portion of Figure [21 smoothed and made negative, with potential 
supergranule features shown by the lighter colors. The identified cell centers are 
marked by the black asterisks. Figure [3] (right) shows the cell centers (now de¬ 
noted by the white asterisks) plotted over the same portion of the corresponding 
unsmoothed travel-time map. The results look reasonable. Then, a cutout is 
taken centered at each identified position of radius ~ 50 Mm and stacked and 
averaged. This is done for each travel-time map for each filter. It is then carried 
out for each of the inverted velocity maps using the same coordinates. 

On average, the algorithm identihed 625 cells per date, for a total of ~ 20,000 
over the entire 32 day period. Though larger than th e number (^^ 5, 500 on 
account of selecting only the largest ones) found by ISvan d^ (|2012h . this is 
substantiall^,smalle^^hai^^hejiumbe^_[^_55j,000]_average^ove^^^Duvall and 
Hanasoge (l2013l 'l~ since we are using data rather close to disk center. However, 
a sample of 20,000 cells should be sufficiently large for our purposes and serves 
to reduce the level of noise in the travel-time and flow maps by a factor of 
of cells = 142. 


4. Results 

4.1. Travel-Time Averages 

Figure 0] shows the resulting supergranule feature derived by averaging the ‘oi’ 
travel-time maps about all ~ 20,000 cell-center locations for each of the eight 
filters. Maps from filters tdi”td 4 show a well-defined circularly-symmetric neg¬ 
ative travel time at cell center, surrounded by a positive ring arising from the 
contribution of adjacent supergranules. Looking to filter tds, some interesting 
structure begins to arise - a small positive signal emerges at the center of the 
cell. In fact, all subsequent Liters tde-tdg exhibit this same central positive ‘oi’ 
travel time, which appears to rise in magnitude and peak ^ 1 sec. A large 
(spatially) ± travel-time asymmetry in the east-west direction is observed in 
the ring surrounding the supergranule for phase-speed Liters td 4 -td 8 , though its 
origin is unknown. 

Often, positive and negative ‘oi’ travel times are interpreted as regions of 
inLow and outLow respectively. This assumption is well-warranted in the shal¬ 
low subsurface layers where low phase-speed travel times are more signiLcantly 
effected by horizontal Lows (and where horizontal supergranular Lows are likely 
largest). It is worth noting, however, that waves propagating to deeper layers 
(waves of higher phase-speed) travel along paths that are more vertical when 
leaving and approaching the solar surface, and are therefore likely to be more 
sensitive to the vertical Lows within a cell. Filters tds-tds isolate signal from 
waves which propagate quite deeply through the solar interior, with lower turn¬ 
ing points ranging from 7-20 Mm. It is therefore important not to necessarily 
interpret these deeper positive and negative ‘oi’ measurements strictly as inLows 
and outLows. To see more clearly what is going on at these depths, it is necessary 
to examine the averaged Low maps. 
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Figure 4. HMI ‘oi’ travel-time maps for each of the eight phase-speed filters tdi-tdg after 
averaging over 20, 000 supergranule features. Filter number is shown in the upper left-hand 
corner of each panel with corresponding parameters given in Table [2] 


4.2. Velocity Averages 


Figure [5] shows similar cuts taken through the inverted horizontal flow compo¬ 
nents of the average supergranule. The inversion depth for which each average 
was computed is shown in the upper left-hand corner of each panel. The root- 
mean-square (RMS) average of all non-zero flow elements in each map is shown 
in the upper right-hand corner, with the associated noise level given in the lower 
left-hand corner. The noise e stimates for the first six inversion depths are given in 
Table 3 of IZhao et all (120121 ). and have been corrected here to take into account 
the statistical averaging. As they are not explicitly stated in the table, the last 
two inve rsion depths have been assigned the same noise level as the deepest layer 


given in Zhao et 'al\ i 2012h . 


The first four panels exhibit a nearly symmetric radial outflow down to a 
depth of 5-7 Mm, roughly consistent with some of the earlier studies mentioned 
in Section [1] At a depth of 7-10 Mm, the sign of the flow changes around cell 
center, and an inflow appears. As the remaining panels show, this extremely 
weak inflow is seen to subsist at all subsequent depths. Interestingly, the near¬ 
surface (0-1 Mm, panel 1) RMS flow amplitudes are substantially lower (by 
an order of magnitud e) than the often-quoted values of ^ 250-350 ms~^ for 
super granulation (e.g. iHartl . [l954 : Simon and height^ . 1964 Hathaway et all 
2OO2II . 

Statistically aver aging the flow field enables u s to make direct comparisons 
with model DH 2 of Duvall and Hanasogel ( 2013 ). The maps in Figure [S] were 
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Figure 5. HMI {vx-,Vy) flow maps for each of the eight inversion depths defined in the text 
(shown in the upper left-hand corner of each panel) after averaging over ~ 20, 000 supergranule 
features. The RMS value of all non-zero flow elements in each flow map is shown in the upper 
right-hand corner of each panel, with the associated noise level given in the lower left-hand 
corner. 


stacked at each central depth, and the horizontal divergence (Vh • rih) was com¬ 
puted from the flows. Figure [S] shows a cut in depth through the divergence 
at y = 0 for both the HMI supergranule (top) and DH 2 (bottom). The HMI 
divergence has been interpolated onto the DH 2 grid for comparison. We find 
that the HMI cell exhibits an outflow which peaks at a depth of about 2 Mm, 
and extends to approximately twice the depth of the DH 2 outflow. The HMI 
supergranule shows the transition region at a depth of ^ 8 Mm, with a weak 
extended inflow reaching down to ^ 20 Mm. As expected from the vector map 
RMS values, the HMI divergence is an order of magnitude smaller than that of 
DH2. 

Figure |7| shows one-dimensional depth cuts through the Vx flow component for 
the HMI and DH 2 supergranules. These profiles were taken through the points 
of maximum surface Vx- Figure |S| shows a cut through the HMI supergranule in 
vector form. The contours mark the 20, 40, and 60 ms“^ scalar velocity levels. 
This can be directly compared with DH 2 in Figure |TJ 
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Figure 6 . Cuts in depth at y = 0 through the horizontal divergence (Vh • for average 
HMI supergranule (top) and DH 2 (bottom). 



Figure 7. One-dimensional cuts in depth through the average HMI supergranule and DH 2 
Vx- Also shown are cuts through the best-fit flow models to all short and large-distance ‘oi’ 
travel-time measurements (BFMi and BFM 2 respectively). These models are discussed in 
Section 15.41 All profiles are taken through the points of maximum surface Vx for the flow 
fields. 


4-2.1. Weak Near-Surface Flow Amplitudes 

It could be argued that the small horizontal flow amplitudes observed near the 
surface for the HMI supergranule are simply a result of the statistical averaging 
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Figure 8. A cut in depth at y = 0 through the average HMI supergranule {vx,Vz) flow com¬ 
ponents. The contour lines denote the 20, 40, and 60 ms“^ scalar velocity levels corresponding 
to the color bar. 


procedure. If the supergranule cell centers have not been correctly aligned (say 
they are shifted by a pixel or two from one another), we might expect the overall 
flow amplitude to diminish significantly after combining 20,000 features. To 
check if this was indeed the case, the near-surface horizontal flow magnitude 
(timag = Vy) was computed and averaged in an azimuth about cell center, 

individually for each supergranule in our sample. This gives an idea for what 
typical supergranule flow amplitudes are that is independent of the overall sta¬ 
tistical averaging scheme. Figure IHl shows a histogram of the resulting maximum 
values for Vmagir, z = 0) for every supergranule. We find that the majority of 
HMI supergranules have maximum flow amplitudes in the range of 60-80 ins“^, 
which is indeed significantly less than values quoted elsewhere. 

It is also possible that since our analysis focused only on data within ±24° of 
disk center, we were not truly “capturing” the full strength of the supergranule 
horizontal flow f i eld lik e we would if we were looking close to the solar limb (e.g. 
[Hathaway et ai . 2002h . In addition to the twenty-five data patches discussed in 
Section[2l each date also included another four patches centered at (0, ±68° ) and 
(±68°, 0) in latitude and longitude. Flow magnitudes were examined individually 
for the supergranules identified in these patches over several days and compared 
to the ones near disk center. We find no major quantitative differences in terms 
of overall flow amplitudes for the sample when looking closer to the limb. 

Another factor contributing to the reduced flow amplitudes likely comes from 
the smoothing inherent in the inversion procedure itself, since what is recovered 
is essentially a convolution between the inversion averaging kernel and the flow 
field at some target depth. The resulting flow map therefore has a resolution 
approximately corresponding to the width of the averaging kernel (generally 
considered to be the averaging kernel FWHM in the case of a Gaussian-shaped 
function). In the best-case inversion scenario, the averaging kernel FWHM would 
be as narrow as possible, though this is generally not practical due to the re¬ 
sulting magnification of noise in the problem. As such, this constraint is often 
relaxed so that near-surface inversions ty pically have resolutions < 10 Mm or so. 
IPeGrave. Jackiewicz. and Remnel ( 2014 1 and some further testing showed that 
the convolution of a Gaussian function of FHMW = 10 Mm with a diverging 
flow field with 250 ins“^ flows can reduce the amplitude by about 35%. Figure ITUl 
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Figure 9. Histogram showing the maximum value of the azimuthally-averaged near-surface 
flow magnitude (umag = \/'^i every supergranule in our sample. 


shows cuts through the x, y, and z components of the HMI averaging kernel for 
the near-surface (0-1 Mm target depth) Vx how inversions. The kernel is found 
to peak at a depth of roughly 2 Mm, and is quite broad in depth relative to 
inversions performed elsewhere ('e.g. lJackiewicz. Gizon. and Birchl . l2008l : ISvandal . 
2012 : DeGrave. Jackiewicz. and Remodl . I20l4i . The effective horizontal and 


vertical resolutions of the kernel are roughly 10 Mm and 4 Mm respectively 
(bottom row figures), and so we would expect some influence from this factor. 
The effect was tested by convolving the averaging kernel with the model DH 2 
flow field. This resulted in an RMS Vx flow amplitude that was 50% lower than 
the actual model value at a depth of 0.5 Mm. 

For comparison, additional cuts through the HMI averaging kernels for inver¬ 
sions centered at depths of 1-3, 3-5, and 5-7 Mm are also shown in the bottom 
right panel of Figure ITOl We find that the three near-surface inversion averaging 
kernels are almost identical in terms of their sensitivity in the z-direction, and 
also the depth at which peak sensitivity occurs. This likely explains much of the 
observed correlation between the shallowest flow maps in Figure [5] 

Another likely effect is the use of standard phase-speed filters employed in the 
HMI time-distance pipeline fTableEl) . iDuvall and Hanasogd ( 201^ demonstrated 
that these rather narrow (in phase space) filters do not capture the full scattered 
wavefield from supergranules, causing travel times to be smaller than what would 
otherwise be found if wider filters were used. To this end, they advocate the use 
of very wide {A£ = 400) phase-speed filters. 
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Figure 10. From left to right: Cuts through the y, and 2 components of the HMI averaging 
kernel for the near-surface (0—1 Mm target depth) Vx flow inversions. The square panels in the 
top row show each averaging kernel component after integrating over depth. The middle row 
panels show cuts in depth through the components along y = 0. Red and black contour lines 
denote 50% of the maximum value of the x-component of the kernel, respectively, while the 
blue and green curves mark the ±5% contours respectively. The left-most bottom panel shows 
normalized ID cuts through the depth-integrated x-component of the averaging kernel along 
y = 0 (black curve) and x = 0 (red curve). The right-most bottom panel shows a normalized 
ID cut in depth through the x-component of the averaging kernel along y = x = 0 (black 
curve). Also plotted here for comparison are ID cuts in depth through the x-component of the 
HMI averaging kernels for inversions at depths of 1-3, 3-5, and 5—7 Mm (colored curves). 


5. Comparisons Through Forward Modeling 


The main observational result of the Duvall and Hanasoge ( 20131) paper was the 
fact that measured ‘oi’ travel times at the center of an average supergranule 
over the distance range A = 2-24° approached a nearly-constant 5.1 sec for 
the largest A values. Travel times calculated for a series of models in the ray 
approximation were compared to these measured travel times, where it was 
determined that model DH 2 provided the best fit to the data. Later work by 
Duvall. Hanasoge. and Chakrabortvl ( 2014 ) also compared large-distance mea¬ 
sured and ray-approximation travel times computed in the ‘we’ and ‘ns’ quadrant 
configurations, again concluding that DH 2 adequately represented the data. 
The application of ray theory to their work was validated through favorable 
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Figure 11. Comparing the statistically-averaged HMI supergranule ‘oi’ travel times (HMI 
Avg. SG oi) to forward-modeled Born-approximation travel times for the average HMI su¬ 
pergranule (HMI Avg. SG Born) and model DH 2 . These forward-modeled travel times were 
computed by convolving each supergranule feature with the HMI time-distance pipeline Born 
kernels. All ‘oi’ travel-time values were taken from cell center. 


comparisons of ray-theory travel-time calculations and travel times measured 
using simulated data having a prescribed vertical flow field. 

In this section, we further test this flow model through forward modeling in the 
Born approximation. Forward-modeled travel times are compared to measured 
ones using kernels produced using two very different sets of phase-speed hlters, 
each of which is valid over a different, yet overlapping range of distances. The first 
is a short-distance range spanning A = 0.5-4.5°, and will be collectively referred 
to as Ai hereafter. The second, a large-distance range spanning A = 2-24°, will 
be referred to as A 2 . 

5.1. Short-Distance Born Comparison Using Standard HMI Phase-Speed 
Kernels 

To begin, the actual set of Born kernels used in the HMI time-distance pipeline 
was downloaded from JSOC. These kernels were produced using the set of stan¬ 
dard HMI phase-speed hlters whose parameters are given in Table [U and are 
valid over Ai. These were convolved separately with the average HMI supergran¬ 
ule and DH 2 how helds to produce a series of forward-modeled travel-time maps. 
The central ‘oi’ travel-time value from each map was plotted as a function of A 
for both features, as shown in Figure [TTJ Also plotted here for comparison are 
the measured travel-time values of the statistically-averaged HMI supergranule 
(those shown in Figure S]). 

It is apparent that the forward-modeled travel times of the average HMI 
supergranule are somewhat inconsistent with the corresponding statistically- 
averaged measurements over distances A < 2°, differing in amplitude by a factor 
of 3.5 at the smallest A values. At larger distances, agreement is closer, though 
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Figure 12. Left: Comparing our large-distance ‘oi’ (left) and quadrant (right) Born-approxi¬ 
matio n forward-modeled travel times (dashed lines) to the corresponding iDuvall and Hanasogd 
ll201.'jl measured (D&H curves) and ray-approximation times (solid lines) for model DH 2 . 


the measured travel times are small and positive, while the forward-modeled ones 
remain negative. At short distances, DH 2 exhibits travel-time amplitudes that 
are larger (more negative) than those of the statistically-averaged measurements 
by a factor of 1.5-2. The stati stically-averaged measuremen ts in this figure com¬ 
pare favorably to Figure 3 of iDuvall and Hanasogel ( 2013h . where travel times 
for an average supergranule were measured from HMI data after first applying 
a series of nominal filters similar to those used in the HMI pipeline. 

The comparison between the measured HMI and forward-modeled travel times 
using the average supergranule flow field is essentially only a test of the inver¬ 
sion algorithm. To the extent that some of the inputs to the inversion, such as 
the travel-time noise covariance matrices and some of the constraint parameter 
values, are unknown, the agreement is not unreasonable. 


5.2. Large-Distance Born Comparison Using Wide Phase-Speed Kernels 


As a secondtesLa^^et_of^orma£2£2SSi2ii2IL2£IiSiilii2j^££ii£l§J2LJS22J®i^'‘^^ 


and Gizon, 200?! ) was comp uted over A 2 using the sam e fourteen wide phase- 


speed filters empl oyed in the Duvall and Hanasogd (201 3ll and Duvall. Hanasoge, 
and Chakraborty ( 2014h studies. Such large-distance, wide-filter Born kernels 
have not previously been used before in a time-distance analysis. These kernels 
were used with the DH 2 model to provide forward-modeled travel-time maps. 
Figure 1121 compares the ‘oi’ ( left) and quadrant (right) cell-center values to the 
Duvall and Hanasoge ( 2013h and iDuvall. Hanasosre. and Chakraborty ( 20141 ) 
measurements and ray-approximation times. Though they differ in amplitude, 
below A = 3.5° the Born and ray-approximation ‘oi’ travel times show similar 
trends, initially rising sharply and peaking around A = 2.5° before falling off. In 
the range A = 7-8°, something interesting happens: rather than approaching a 
value ~ 5.1 sec, the Born travel times exhibit a second large peak (of unknown 
origin, as visual inspection of the kernels shows no obvious anomalies) before 
decreasing again and leveling off at values ~ 2 sec at larger distances. Unlike 
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Figure 13. Top row: Horizontal and vertical contributions to all large-distance DH 2 Born-ap¬ 
proximation forward-modeled ‘oi’ (left) and quadrant (right) travel times shown in Figure [12] 
Bottom row: Cuts at y = 0 in the east-west direction through the A = 11.76° DH 2 ‘oi’ 
(left) and ‘we’ (right) Born-approximation travel-time maps with the individual horizontal 
an d vertical contributions to these time s inclu ded. These can be directly compared to Figure 2 
of iDuvall. Hanasoge. and Chakrabortvl ll20l4l . The starred values in the bottom row figures 
correspond to the starred points in the top row figures. 


the ray travel times, it is apparent that these Born travel times do not closely 
match the data at any distance. Examination of the quadrant travel times also 
shows large discrepancies. Here, we find that the Born travel times are small and 
relatively constant over all distances, differing from both the measurements and 
ray-approximation times by an order of magnitude at small A, and a factor of 
4.5 at larger distances. 

5.3. Disagreement Between Measured and Forward-Modeled Travel Times 

Based on these comparisons, it is difficult to draw a firm conclusion about 
whether or not DH 2 is indeed a reasonable representation of an average su¬ 
pergranule. It could be that the model is inadequate in light of the disagreement 
between measured and Born travel times over these distance regimes, or it could 
be that the Born kernels used in this analysis (those obtained directly from the 
pipeline, as well as the newly-computed large-distance ones) are not correct. 

In terms of the wide-filter results over A 2 , some disagreement between the 
measurements and our forward-modeled travel times may be expected, as the 
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computation of these kernels assumes a plane-parallel geometry, neglecting the 
curvature of the solar surface - an assumption which admittedly begins to break 
down at large distances. It is surprising, however, that the Born travel times do 
not match the data well even at relatively short distances where the deviation 
from plane-parallel geometry is still quite small. 

It is also possible that some mismatch is due to the inability of these kernels 
to adequately separate the horizontal and vertical flow components at these 
distances, introducing a type of ‘cross talk.’ To test this, we computed the 
horizontal and vertical contributions to the ‘oi’ and ‘we’ forward-modeled travel 
times for each distance. These are shown in the top left and top right panels of 
Figure [T3] respectively. These figures show that the kernels are in fact able to 
separate the contributions fairly well - in other words, one component contributes 
to the signal significantly more than the other - particularly for the ‘oi’ times, 
except in the range A = T-S"”. The bottom left and bottom right panels of 
Figure [13] show east-west cuts at ?/ = 0 through the ‘oi’ and ‘we’ maps for 
one particular distance (A = 11.76°), along with the individual horizontal and 
vertical contributions to these travel times. The starred points in the bottom 
row hgures correspond to the starre d points in the top row hgures. We hnd that 
the ray kernels used in Figure 2 of iDuvall. Hanasoge. and Chakrabortvl ( 2014l l 
separate the contributions more fully at these distances than these particular 
Born kernels do. 

The computation of the Bor n kernels relies on an accurate mode ling of the 
HMI acoustic power spectrum ( Birch. Kosovichev. and Duvalj BoOdh . This can 
be quite challenging in some cases, particularly when usin g filters of higher phase- 
speed, and the wide filters of iDuvall and Hanasogel (l2013h . It is therefore possible 
that some discrepancy, at least over A 2 , is coming from an insufficient modeling 
of this data power. To get a sense for how well the power has been modeled in each 
case, we compared the filtered data and model power spectra after integrating 
separately along frequency and wavenumber. Figure 1141 shows the integrated 
data power (solid lines) and model power (dashed lines) plotted together for 
each kernel. Filter number is given in the upper right-hand corner of each panel. 
The power has been normalized to unity in each case for easier comparison. 
We note that the same phase-speed filters are used for both the model and 
data power spectra. To explore this idea a bit further, we examined our set of 
wide-filter kernels more closely. 
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Figure 14. Filtered data (solid lines) and model (dashed lines) power spectra after integrating 
over frequency (a) and wavenumber (b) for all large-distance Born kernels used in this analysis. 
The power has been normalized to unity in each case for easier comparison. Filter number 
is shown i n the upper right-hand corn er of each panel. These are the same fourteen filters 
used in the lDuvall and Hanaso^ (1201311 study, with larger filter numbers correspond to larger 
phase-speed values. They correspond to wave travel distances of about 2—24°, in approximately 
2 ° bins. 
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In terms of power integrated over frequency ('FigurefUk'). the model power is 
found to deviate most significantly from the data at high wavenumber for the 
lowest phase-speed filter, tdi. As we look to filters of higher phase-speed, the 
model power appears to quickly “pull away” from the data, with the fit becoming 
progressively worse at nearly all wavenumber values. Beyond filter td 4 , the model 
profiles generally are indistinguishable from one another. In most cases, we are 
unable to model the power at low wavenumber when using these wide filters. 

To some extent, the agreement between data and model power appears to be 
slightly better when integrated over wavenumber (Figure [HJo), though we do 
have difficulty getting good agreement between the spectra at higher frequency 
in most cases. We find that the profiles for each of the filters exhibit a strong 
oscillatory behavior, which is also present in Figure [Hk. though to a somewhat 
lesser degree. Increasing the number of frequency and wavenumber model grid 
points usually helps to counteract this, however this can slow kernel computation 
time down considerably, and, after a certain point, starts to become impracti¬ 
cal. We have currently not tested what overall effect this behavior has on the 
computation of forward-modeled travel times. 

These types of figures are useful for understanding modeling uncertainties and 
are not always provided in publications. Unfortunately, we do not have available 
the HMI pipeline model power spectra with which to compare. 


5.4. The Best-Fit Flow Models 


In light of the fact that the best iDuvall and Hanasogd (1201311 flow model is in 
disagreement with the forward travel times in the Born approximation, we were 
interested in seeing what kind of supergranule features could provide the best 
fit to the measurements made over the separate Ai and A 2 regimes, under the 
assumption that both sets of Born kernels are accurate . Usin g the continuity 
equation and a Model S ( Christensen-Dalsgaard et al . 1996h density profile. 


512, 000 Duvall & Hanasoge-style flow models (i.e. mass-conserving models with 
a Gaussian flow profile) of varying peak amplitude, peak Vz depth, and 
Vz FWHM were produced in a Monte Carlo-like fashion. Imposing observational 
constrains on these models to have a surface Vz{z = 0) ~ 10 ms“^ and a surface 
= 0) = 100-500 ms“^ resulted in a total of 590 “viable” flow models. 

Forward-modeled Born travel times were computed for these flow models over 
Al and A 2 . The best-fit models to all short and large-distance ‘oi’ measurements 
(in a least-squares sense, referred to as BFMi and BFM 2 , respectively, hereafter) 
were identified, and have the following parameters: 


BFMi: 

Vz{z = 0) = 9 ms-i and vf^^{z = -1.7 Mm) = 154 ms-\ FWHM^ = 1.7 Mm 
v-h{z = 0) = 148 ms“^ and v^^^{z = —1.2 Mm) = 559 ms“^ 

BFM2: 

Vziz = 0) = 14 ms-i and u““(z = -3.6 Mm) = 352 ms-\ FWHM^ = 3.3 Mm 
Vhlz = 0) = 135 ms“^ and v^^^{z = —2.5 Mm) = 786 ms“^ 


SOLA: Paper_solar.tex; 8 April 2015; 0:26; p. 19 













K. DeGrave and J. Jackiewicz 




Figure 15. Short (left) and large-distance (right) ‘oi’ travel-time comparison for all flow 
features. These plots are identical to Figures 1111 and 1121 only they now include the BFMs 
(dashed-star lines), RSGi, and RSG 2 (these features are discussed in Section [6jl values for 
comparison. The measured travel times in the left and right-hand figures are denoted by the 
‘HMI Avg. SG oi’ and ‘D&;H (2012) oi’ curves respectively. Recall the left (right) plot implies 
the use of narrow (wide) phase-speed Alters. 


Depth cuts through the two models are shown in Figure [7] along with DH 2 
and the average HMI supergranule. BFMi is shallower and of lower amplitude 
than DH 2 , while BFM 2 is of larger amplitude and is much more extended in 
depth. 

The resulting forward travel times from these models are shown in Figure [15] 
These plots are identical to Figures [TT] and [TTl only they now have the additional 
BFM values included for comparison (dashed-star lines). We find that the models 
do in fact provide better fits to the measurements when compared to the DH 2 
Born-approximation travel times, especially with regard to the Ai measurements 
(left). The fit over A 2 (right), however, is still quite poor. Rather than comparing 
only cell-center ‘oi’ travel time values, we also compared the full measured and 
modeled travel-time maps when available. This was only possible over Ai, as we 
did not have available the full Duvall & Hanasoge travel-time maps with which 
to compare. The travel-time maps from each model were compared separately to 
those of the averaged HMI measurements for each filter. For each comparison, 
the relative error was computed by finding the RMS difference between maps. 
Overall, we find no difference in terms of which model best fits the data when 
comparing the full travel-time maps versus using only the cell-center values as 
shown in Figure fT^ (left 1. 

Experimenting with many different flow models like these provides a few inter¬ 
esting conclusions. We find that the model that best approximates the measured 
travel times over Ai does not fit the A 2 , wide phase-speed measured times well, 
and vice versa. It is very difficult to simultaneously get a reasonable fit over both 
distance regimes with a single model, which is not surprising given the different 
filtering. It also appears that the parameter space is somewhat degenerate in the 
sense that there are many flow models for which the fit to the measured travel 
times is nearly equally good. These models are often quite dissimilar in peak Vz 
amplitude, peak depth, and FWHM^. To illustrate this, Figure Uni shows the sum 
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Figure 16. Figures (a), (b), and (c) show the sum of the residuals between the A 2 measured 
and forward-modeled ‘oi’ travel times for each of the 590 flow models as a function of the 
model parameters. These have been normalized to the largest residual value in each panel. 
Small residuals denote flow models that more closely approximate the measurements. These 
figures serve to illustrate the degeneracies between many of the models. Panel (d) shows the 
peak Vx. velocity for each of the models, with values corresponding to the color bar. The color 
bar values have been divided by a factor of 100 for easier viewing. White areas in each panel 
represent regions where no viable models (i.e. ones obeying the surface flow constraints) were 
identified. The circle and star points mark the locations of DH 2 and BFM 2 , respectively, in 
this parameter space. 


of the residuals between measured and forward-modeled ‘oi’ travel times over A 2 
for each of the 590 flow models as a function of these parameters. These have 
been normalized to the largest residual value in each panel, with the smallest 
residual values representing models that are able to more closely match the 
measurements. White areas represent regions where no viable models (i.e. ones 
obeying the surface flow constraints) were identified. The circle and star points 
mark the locations of DH 2 and BFM 2 , respectively, in this parameter space. The 
overall difference between them is not substantial. 


6. Simulated Supergranule-Sized Flow Feature 


It is possible to derive an alternative model of supergranulation from the large- 
scale flow struct ures present in recent realistic magnetohydrodyamic quiet-Sun 
simulations fe.g. iRem uell. 2 014ll. Having already analyzed on e such simulation 
previously (QS2 in DeGrave. Jackiewicz. and Rempell . l2014[) . these data were 
already available. The near-surface flow field of this simulation was smoothed 
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Figure 17. One-dimensional cuts in depth through the RSGi and RSG 2 Vx. (left) and Uz 
(right) flow components. These profiles are taken through the points of maximum surface Vx 
and Vz, respectively, for the two flow fields. 


with a two-dimensional Gaussian function of cr = 4 Mm to help remove small- 
scale fluctuations and accentuate the regions of large-scale outflow. Points of 
maximum divergence were then located from the flows, thereby marking the 
positions of outflow centers. Averages were taken about seven of these features 
to produce a single average “supergranule” (referred to as RSGi hereafter). Guts 
through the Vx and how components of RSGi are shown in Figure flTl (solid 
lines). These can be directly compared to the HMI, DH 2 , and BFM prohles in 
Figure[7l The RSGi feature is clearly very different from the others in its overall 
flow structure, particularly its deeper extent in depth to at least 15 Mm. 

The supergranule obtained this way does not adequately satisfy the continuity 
equation (see Table O RSGi). To determine if this is one of the reasons for the 
different flow profile, we used RSGi as a starting po i nt to derive a fully mass- 
conserving flow model like in iDuvall and Hanasogd ( 2013 ). A model (referred 
to as RSG 2 hereafter) was chosen to give a satisfactory match to the RSGi Vz 
profile in depth, with the Vx profile being determined completely by the choice 
of Vz- Figure [T71 shows cuts through the Vx and Vz flow components of RSG 2 
(dashed lines). The model matches reasonably well the overall amplitude of the 
RSGi flow components as they vary in depth, though does not have the large- 
amplitude spike in Vz that RSGi exhibits near the surface. It also has stronger 
near-surface horizontal flows. 

Forward-modeled travel times were computed for RSGi and RSG 2 in the Born 
approximation over both Ai and A 2 distance ranges, with the results plotted 
in Figure [T51 Over Ai, the two features show similar values at all distances. 
However, they are generally unable to approximate the HMI measurements well, 
being larger in amplitude by a factor of 1.5-2.5 at small A. This is not too 
surprising given the very strong near-surface horizontal flows. We note that 
unlike the statistically-averaged HMI measurements, none of the models produce 
positive ‘oi’ cell-center travel-time differences over the larger distances of Ai. 

Over A 2 (Figure [T51 right), RSGi and RSG 2 travel times show larger dis¬ 
agreement between each other, however, RSG 2 matches the measured travel 
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times at nearly all distances. Like BFM 2 , RSGi and RSG 2 tend to overshoot the 
measurements at shorter distances, and also exhibit the anomalous peak between 
A = 7-8° that the other models also show. Though not shown in Figure [151 the 
‘we’ and ‘ns’ quadrant forward-modeled times for RSGi, RSG 2 and BFM 2 show 
no significant improvement over the DH 2 Born values in Figure |T2j 


7. Conclusions 


HMI time-distance pipeline velocity data products suggest that the “average” 
supergranule exhibits the following structure: 

• A radial outflow down to a depth of 7-10 Mm 

• An transition from outflow to inflow at a depth of 7-10 Mm 

• A weak inflow flow extending down to a depth of ~ 20 Mm 


These proper ties differ significantly f rom th e shallow model of supergranulation 
proposed by Duvall and Hanasoed (2013h. and further supported bv Duvall. 
Hanasoge, and Ghakraborty ( 2014ll . The near-surface RMS velocity of the av¬ 
erage HMI feature is also found to be considerably smaller (by an order of 
magnitude) than the values typically associated with solar supergranulation. 
Much of the overall HMI supergranule structure and its weak flow amplitudes 
can be explained by examining the HMI pipeline averaging kernels for the 
near-surface inversions, which are found to be very broad in depth, and nearly 
identical in terms of sensitivity along the z-direction. Additional reduction in 
flow amplitude could be attributed, at least in part, to the use of narrow stan¬ 
dard phase-speed filters in the analysis, along with the smoothing of recovered 
flows inherent in time-distance inversions. A convolution between the HMI near¬ 
surface Vx inversion averaging kernel with model DH 2 resulted in a 50% reduction 
in Vx RMS flow amplitude. 

To identify possible sources of discrepancy between the data and model, we 
studied forward-modeled travel-time differences in the Born approximation to 
compare with the HMI measurements. Using kernels computed with the standard 
(narrow) phase- speed filters in the HMI pipe line over distances Ai = 0.5-4.5°, 
the model from Duvall and Hanasogel ( 20I3ll . which was originally derived in 
a ray theory approach, does not match the measured HMI travel times for an 
average supergranule fFigure ITT]) . 

At larger distances (A 2 = 2-24°), travel-time measurements exist only from 
Duvall and Hanasogel ()2013l l who employed wider phase-speed filters. Born ker¬ 


nels computed using these Liters also do not reproduce what ray theory predicts 
for the shallow supergranule model (Figure [T^). It is unclear if it is appropriate 
to compute Born kernels with such filters, due to potential modeling inaccuracies 
over such a large region of the power spectrum. For these rather small travel¬ 
time differences, corrections due to the curvature of the Sun that are not taken 
into account here could become important too. 

When we instead assume that the Born kernels computed with short-distance, 
standard phase-speed filters and the large-distance, wide phase-speed filters are 
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correct, we determined a best-fit supergranulation model to the measurements 
over each of those distan ce ranges. The resulting sup ergranule is different in 
each case (Figure [5]) • The Duvall and Hanasogel ( 201, “ill model falls somewhere 
between these two solutions in terms of amplitude and depth variations. 

Another representation of a supergranule was derived by averaging over large- 
scale flow features in a realistic 3D numerical simulation. Minor modifications 
to the flow profile to make it fully satisfy the continuity equation were made. 
The str ucture has very strong flows that extend to much larger depths than the 
HMI or Duvall and Hanasosrel ( 201,3h supergranules. In the Born approximation, 
travel times computed from this supergranule only match the data reasonably 
well over a small range at large distances. _ 

It is quite possible that neither the lDuvall and Hanasoee ( 2013h flow model 
nor the average HMI feature are accurate representations of an actual s upergran¬ 
ule. T he latter may be especially true in light of the recent work by Hanasogg 
( 2014l l where it was found that, even in the case of an ideal inversion scen ario, 
model supergranule flows could not be accurately recovered. Svandal (201 ,^ also 
suggested that inaccuracies in the forward problem may be amplified in the 
inversions process, leading to erroneous subsurface flows recovered through time- 
distance helioseismology. In any case, we And that the supergranules implied 
from three different approaches are not consistent with each other. The one 
most closely in agreement with past surface measurements (by construction), as 
well as large-distan ce travel-time differences inte rpreted from ray theory, is the 
model proposed by Duvall and Hanasogel ( 2013h . More work needs to be done 
to show that the Born approximation can be used to determine this structure. 
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Appendix 

A. Mass Conservation of the Various Flow Features 


The average HMI and RSGi features were checked to see how well their flows 
conserve mass relative to the various models described throughout this work. 
We do this by expanding the continuity equation (term [A] in Equation [T|) into 
its separate horizontal and vertical components (terms [B] and [C] respectively) 
assu ming a time-independent Model S density profile (' Christensen-Dalsgaard 
et ai, 1996l l which varies only in the ^-direction: 


V • ipv) = Vh • (puh) + = 0. 



( 1 ) 
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Table 3. Mass conservation ratios for the 
HMI average supergranule, DH 2 , RSGi, 
RSG 2 , and the BFM. Parameters B and C 
are the spatially-integrated values computed 
from the expanded continuity equation in 
Equation ^ 


Feature 

(B -f C)/B 

(B-hC)/G 

HMI 

-3.55 X 10-1 

2.62 X 10-1 

DH2 

1.54 X lO-'i 

-1.54 X 10-1 

BFMi 

1.36 X lO-'i 

-1.36 X 10-1 

BFM2 

-1.55 X 10-1 

1.55 X 10-1 

RSGi 

-3.04 X 10-2 

2.95 X 10-2 

RSG2 

1.12 X 10-® 

-1.12 X 10-® 


The model supergranules are constructed analytically to be mass conserving, 
though evaluate numerically to some small non-zero value due to numerical 
gradients and machine precision. As such, we check mass conservation by calcu¬ 
lating terms [B] and [C] within a 30 Mm radius of cell center, and integrating 
them along all three spatial dimensions. Their ratios are then taken with respect 
to [B + C]. In the case of a mass-conserving flow, the ratios (B -|- C)/B and 
(B -I- C)/C should be of the same magnitude and of opposite sign. The mass 
conservation ratios for each feature are presented in Table [31 As expected, DH 2 , 
RSG 2 , and the BFM all satisfy continuity to a high precision. The HMI and 
RSGi supergranules, to varying degrees, do not. 
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